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A general covariant extension of Einstein 's field equations 
is considered with a view to Numerical Relativity applica- 
tions. The basic variables are taken to be the metric tensor 
and an additional four- vector Einstein's solutions are re- 
covered when the additional four-vector vanishes, so that the 
energy and momentum constraints amount to the covariant 
algebraic condition = 0. The extended field equations can 
be supplemented by suitable coordinate conditions in order 
to provide symmetric hyperbolic evolution systems: this is 
actually the case for either harmonic coordinates or normal 
coordinates with harmonic slicing. 



I. INTRODUCTION 

General covariance is a key feature of General Rela- 
tivity. At a first look, Einstein's field equations can be 
understood as a set of ten second order partial differential 
equations on the ten unknown metric coefficients g^i, 
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TT (T^^ - - T .g^^) 



However, in General Relativity, like in many other field 
theories, physical solutions are given by an equivalence 
class of field values, related to one another by gauge 
transformations. General covariance means that the 
gauge group in Einstein's theory is that of the general 
(smooth) coordinate transformations 



(2) 



Accordingly, the field equations (1) do not (even may 
not) provide enough information to determine the values 
of the ten unknown coefficients 17^,^. On the other hand, 
in numerical applications we must deal with specific met- 
ric components. It follows that a numerical evolution 
system must include a specification of the coordinates as 
an extra ingredient in order to determine the four kine- 
matical degrees of freedom. Fixing four of the ten metric 
coefficients we choose one specific expression for 17^^ out 
of the equivalence class representing the same physical 
solution. A general covariant evolution system would be 
incomplete and conversely, a complete evolution system 
can not preserve general covariance. 

To be more specific, we will carefully distinguish evolu- 
tion systems, associated with a particular gauge choice, 
from evolution formalisms. The later can be defined as 
a set of equations that apply to an equivalence class of 
solutions, prior to a complete gauge specification. Fol- 
lowing this way, we will be able to draw a line between 



general covariant formalisms, where the equivalence class 
is defined by the full group of coordinate transforma- 
tions (2) and partially covariant formalisms, where the 
equivalence class is defined by any restricted subset of 
coordinate transformations. This means also for instance 
that different evolution systems can be obtained from the 
same evolution formalism, as we will see in the following 
paragraphs. 



A. General Covariant Formalisms 

A good example of a general covariant formalism is 
provided by recent numerical relativity works [1,2] based 
upon a well known classical approach. The evolution 
formalism is given by the original field equations (1), al- 
though the DeDonder [3,4] expression of the Ricci tensor 
is used to write down the principal part, namely 



(3) 



(1) where the box symbol stands for the d'Alembert oper- 
ator on functions and we have noted F'' = g^'^T^^ as 



per 

usual. General covariance is not lost in passing from (1) 
to (3), because we only reordered the partial derivatives 
in the Ricci tensor. 

Is is obvious from a comparison between (3) and the 
wave equation for g^^ that we can obtain a symmetric 
hyperbolic system imposing the well known harmonic co- 
ordinate conditions [3,4]: 



-F^ = 0. 



(4) 



Although the first proofs of well-posedness of the result- 
ing evolution system were well known [5,6], the corre- 
sponding proof for the initial-boundary problem, which 
is highly relevant for Numerical Relativity applications, 
has been given recently [1,7,8]. Different evolution sys- 
tems can be obtained from the general covariant formal- 
ism (3) by modifying the harmonic coordinate conditions 
(4) [9,10]. One can even add arbitrary "gauge source" 
terms to the right-hand-side of (4) to obtain a wide class 
of generalized harmonic evolution systems [11]. 

A different evolution formalism can be obtained by us- 
ing the well known 3+1 decomposition [12,13], where one 
considers the space-time sliced hy t = constant hypersur- 
faces. The line element can be written as 



dt^ + 7y [dx' + (3' dt) {dx^ + (3^ dt), (5) 



where the lapse a and the shift /3* represent the kinemat- 
ical degrees of freedom. The field equations (1) can then 
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be translated in terms of the three-dimensional geometry nate transformations (2) that preserve the time slicing, 

of the slices, namely namely 



{di - Cfj}-ftj = -2 a Kij (6) 

(dt - Cp)K,j = -V,aj + a [^^^Rij - 2Kfj + trK Kij] (7) 

^^'>R-tr{K^) + {trKf =Q (8) 

Vk{KK^-8\trK)={) (9) 

where we have restricted ourselves to the vacuum case 
for simplicity. 

Let us notice that the coordinate gauge freedom is not 
limited in any way by translating the four-dimensional 
(4D) field equations (1) into the 3+1 version (6-9). The 
lapse a and the shift /3* can take arbitrary values, so that 
the four gauge degrees of freedom are still at our disposal. 
Although general covariance is not manifest in the 3+1 
equations (6-9), their solution space is still invariant un- 
der general coordinate transformations (2), because it is 
equivalent to the corresponding sohition space of the 4D 
equations (1). Conversely, in the 3+1 version (6-9) it is 
manifestly clear that one can evolve the dynamical de- 
grees of freedom 7ij(t, x^) from any consistent set of ini- 
tial data {"fij (0, x'^), Kij (0, x^)} using the evolution equa- 
tions (6,7). The remaining ones (8,9) can be interpreted 
as constraints. This diversity in the evolution proper- 
ties of both sets of equations is not obvious in the 4D 
version. To summarize: both the 4D equations (1) and 
the 3+1 equations provide equivalent general-covariant 
formalisms, although general covariance is apparent only 
in the 4D version, whereas the evolution properties are 
manifest only in the 3+1 version (6-9). 

B. Extending solution space: taking constraints out 

The 3+1 formalism (6-9) is specially suited for Numer- 
ical Relativity applications. In this context, one usually 
takes advantage of the fact that the energy and momen- 
tum constraints (8,9) are first integrals of (6,7). This 
allows us to enforce the constraints (8,9) on the initial 
and boundary data only, or even to use them to monitor 
the accuracy of the time evolution. But the constraints 
are not enforced by the time evolution algorithm for in- 
terior points, which in its simplest form is based only 
on (6,7) (free evolution approach [14,15]). This "uncon- 
strained" evolution formalism, although perfectly con- 
sistent [16,17], does introduce a strong discrimination 
between the two sets of equations (6,7) and (8,9) that 
breaks the general covariance of the 3+1 formalism. 

To verify this, let us notice first that by replacing 
the full 3+1 formalism (6-9) by the subset (6,7) we are 
actually extending the solution space so as to include 
constraint- violating pairs {7^,!^^}. Now, as the re- 
stricted set of equations (6,7) corresponds to the space 
components of (1), the extended space of solutions will 
be invariant only under the restricted subset of coordi- 



t' = h{t) (10) 
y' = f{x\t). 

This confirms that general covariance is broken in the 
unconstrained evolution formalisms. 

Let us remark that the extension of the sohition 
space of (6-9) is a rule, not an exception, among the 
new formalisms arising after the seminal 1983 work of 
Y. Choquet-Bruhat and T. Ruggeri [18], which opened 
the door to the use of arbitrary shift choices in hyper- 
bolic evolution systems. The bottom line is that the 
constraints (8,9) contained in the original 3+1 evolution 
formalism are at odds with hyperbolicity to the intent 
that some extension is needed in order to modify the 
mathematical structure of the formalism without loos- 
ing the physical solutions. The easiest way of doing this 
is just taking the constraints (8,9) out of the system. 
This is the basic ingredient, although this crucial point 
can be masked by other manipulations on the evolution 
equations, like taking an extra time derivative [18,19], or 
an extra space derivative [20], or using the constraints 
to modify the evolution equations (6,7) [21 26]. The 
resulting formalisms, when supplemented with suitable 
coordinate conditions, provide hyperbolic evolution sys- 
tems that can be used in Numerical Relativity applica- 
tions. From our point of view, these formalisms can also 
be interpreted as providing many non-equivalent ways of 
extending the solution space of (6-9) with at least two 
related common features: constraint equations (8,9) are 
left out of the final evolution formalism and general co- 
variance is broken as a result, even before a specific co- 
ordinate system is selected. 

C. Extending solution space: extra dynamical fields 

A completely different way of extending the sohition 
space is to introduce extra dynamical fields, independent 
of the metric and its derivatives, into the evolution for- 
malism. This alternative has been independently used 
by many Numerical Relativity groups in different ways 
[27-30]. The key idea in these works was to introduce 
three supplementary dynamical fields whose evolution 
equations were obtained by using the momentum con- 
straint (9) . As far as these works were focused on Numer- 
ical Relativity applications [31-33], the supplementary 
quantities were introduced in an "ad hoc" way, break- 
ing even the 3+1 covariance (10) of the formalism. Only 
very recently [34] the same idea has been implemented 
in a way which is at least invariant under the restricted 
subset of coordinate trasformations (10): the extra quan- 
tities are given by a three-dimensional "zero" vector Zj 
which vanishes for Einstein's solutions. During numer- 
ical evolution, however, non-zero values of Zi arise due 
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to truncation errors and the resulting numerical codes 
actually deal with an extended set of solutions. 

Even with this improvement, general covariancc is still 
broken for two different reasons. First of all, although 
the momentum constraint (9) has been incorporated into 
the formalism as the right-hand-side of the time evolution 
equation for Z,, the energy constraint (8) is still taken out 
of the time evolution algorithm. It is obvious that the ex- 
tension of solution space to " energy constraint violating" 
modes can not be invariant under the general coordinate 
transformations (2). This reason alone could be easily 
overcome by proceeding along the lines sketched in [36], 
where every constraint is incorporated into the system by 
adding an extra "Lagrange multiplier" quantity: every 
extra quantity could then be coupled with the previous 
ones in many different ways so that many more new ar- 
bitrary parameters will appear in the resulting evolution 
formalism. But this would just reinforce the second cause 
of general covariance breaking: the lack of covariance of 
the final set of extra quantities. If one wants to extend 
solution space by adding extra fields without breaking 
general covariance, then this set of extra fields should be 
equivalent to some set of well defined space-time quanti- 
ties, independent of the time slicing considered. 



II. GENERAL COVARIANT EXTENDED 
EVOLUTION FORMALISMS 

A. The extended field equations 

We propose to extend the field equations (1) in a gen- 
eral covariant way by introducing an extra four-vector 
Z^, so that the set of basic fields will consist into the 
pair {Qfjiu, Zfj}. The original field equations (1) will then 
be replaced by 

R^u + V^Z, + V.Z^ = 8 TT (V - i T g^,). (11) 

The solutions of the Einstein's solutions can be easily 
recognized among the extended set as those satisfying 
condition [35] 



(12) 



so that the four-vector will provide a simple way to 
monitor the quality of numerical simulations or any other 
kind of approximation scheme. Notice that equations 
(11) are of mixed order: second order in the metric com- 
ponents g^„, but only first order in the extra vector field 
Z^i . This means in particular that terms containing first 
derivatives of Z^ belong to the principal part and that 
they are then relevant to the causal structure of the re- 
sulting evolution systems, as we will see later. 

In order to fully understand the evolution properties 
of the extended equations, let us translate the manifestly 



covariant form (11) into the 3+1 language (5). The co- 
variant four-vector will then be decomposed into its 
space components Zi and the normal component 
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(13) 



where is the unit normal to the t = constant slices. 
The 4D equations (11) can then be written in the equiv- 
alent form: 



{dt - CpYdj = -2 a Ki. 
{dt - C0)Kij = -Viaj 4 



- 2 Kfj + (trK - 2 0) K,^] 



(14) 

[^^'>R^J+V^ZJ+VjZi (15) 

(16) 



{dt - £/3)e = - f^R + {trK - 2 9) trK 

- tr{K^) + 2 VfcZ'^ - 2 {ak/a) Z''] 

{dt - C0)Z, = a [Vfc {K,'' - 6iHrK) + d^Q (17) 

- {ai/a) e - 2 if> Zk] 

where stands for dia and we have restricted ourselves 
to the vacuum case again. 

The evolution properties of the formalism are transpar- 
ent in the 3+1 version (14-17). Only evolution equations 
appear there, without constraints. One can not omit 
any of the equations from (14-17) because all them are 
needed to evolve the extended set of ten dynamical fields 
{7ij, Q, Zi}. This means that general covariance will not 
be broken in Numerical Relativity applications because 
all the equations (14-17) must be used on equal footing 
in the main algorithm evolving the dynamical fields in 
time: there is no room left for equation discrimination. 



B. Recovering Einsteins solutions 

The algebraic condition (12) is useful to check a pos- 
teriori whether a given solution of the extended system 
(11) is actually a solution of the Einstein's field equa- 
tions (1). But it is interesting as well to know a priori 
the necessary and sufficient condition for a given set of 
initial data to generate a physical solution. In this sense, 
one can take the divergence of the extended equations 
(11) to get, allowing for the contracted Bianchi identi- 
ties, 



az^ + R^^ Z"" = 0. 



(18) 



This homogeneous second order equation in Z^ ensures 
that any deviation from the original Einstein equations 
(1) propagates through light cones and also that a suf- 
ficient set of conditions for the initial data to provide 
physical solutions is given by 



Z^{0,x') = dtZ^{0,x') = 0, 



(19) 



where the second equation, allowing for (16, 17), repre- 
sents imposing energy and momentum constraints (8,9) 
on the initial data. 
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This means that the algebraic constraint (12) by it- 
self is not a first integral of the extended system (11). 
Equations (8, 9) appear here as auxiliary conditions so 
that the full set (8, 9, 12) is preserved by time evolution. 
From the practical point of view, this means that one can 
take any set of consistent initial data of Einstein's equa- 
tions (1) and use it with a zero initial value of to get 
an initial data set for the extended equations (11) that 
will generate precisely the same solution. The general- 
covariant gauge-independent equation (18) can be then 
interpreted as a useful tool to understand the propaga- 
tion of the initial constraints (19). 

III. COORDINATE CONDITIONS AND 
SYMMETRIC HYPERBOLIC EVOLUTION 
SYSTEMS 

As stated in the introduction, the evolution system 
is not complete until one provides coordinate conditions 
to fix the four kinematical degrees of freedom. We will 
consider here two different coordinate conditions leading 
to a symmetric hyperbolic evolution system. We will 
start in both cases from one of the two equivalent versions 
of the extended general-covariant formalism (11). 

A. A 4D evolution system in harmonic coordinates 

Let us use the DeDonder [3,4] expression of the Ricci 
tensor to write down the principal part of (11), namely 

-□ 5m. + 5^(r. + 2Z,) + S,(r^ + 2 Z^) = ... (20) 

It is obvious from a comparison with the wave equation 
for Qfjtv that we can obtain a symmetric hyperbolic system 

□ 9t,. = - (21) 

provided that we kill the additional terms in (20) using 
the following extension of the well known harmonic co- 
ordinate conditions: 

ax'' = -V" = 2 . (22) 

The four-vector can be interpreted in this context as 
a sort of " gauge source" , along the lines sketched in rcf. 
[11]. 

Notice that, in contrast with the classical approach, the 
compatibility between the reduced system (21) and the 
coordinate conditions (22) is not an issue in the present 
context because there are now fourteen independent com- 
ponents of the fields {5^^,^^} to be determined by the 
fourteen equations (21, 22). A straightforward analy- 
sis in the 3-1-1 framework shows that this is actually the 
case: the lapse and shift evolution is provided by equa- 
tion (22) and the evolution of the remaining ten degrees 
of freedom, including the four-vector Z^, is given by (21). 
Notice also that equation (22) actually coincides with the 



classical harmonic coordinate condition (4) for physical 
solutions, where vanishes. This is why we talk about 
" harmonic coordinates" to refer also to condition (22) in 
the present context. 

B. A 3-|-l evolution system with harmonic slicing 

Harmonic coordinates are not flexible enough to be 
used in most Numerical Relativity applications. A more 
suitable choice is the "harmonic slicing", in which the 
time coordinate is again assumed to be a harmonic func- 
tion, but the space coordinates are chosen so that the 
mixed components goi vanish (normal coordinates). We 
propose here to keep in the extended case the time com- 
ponent of (22), that is 

goi = , □ a;0 = - 2 ^° (23) 

which can be translated into the 3-1-1 language as 

/3' = , a* Ina = -a {trK - 2 6). (24) 

The principal part of (15-17) can be written as 

dtKij+dk[aX''ij] = ... (25) 

dtZ, + dkHS'^i {trK - 9) - K''i)] = ... (26) 
dtO + dk[a {D'' - E'' - Z'')] = ... (27) 

where we have noted 

Ak = dkO-na), Duj = ^ dk^ij (28) 

A'ii = + ^(A, +D,-2 E, - 2 Z,) (29) 

+ ^{Ai + Di-2Ei-2Zi) (30) 

and Dk = Dkij, Ek = Dijk. 

One can get a fully first order system in the usual way, 
by considering both and Dkij as independent addi- 
tional quantities. Their evolution equations can be then 
easily obtained by differentiating (14) and (24), namely 

dtAk + dkHtrK - 2 6)] = (31) 
dtDkij+dk[aKij]=0 (32) 

The full set of basic independent quantities is then given 
by {a,jij,Kij,Q,Zi,Ak,Dkij} and the non-trivial prin- 
cipal part of the corresponding evolution systems is given 
by equations (25-27, 31,32). 

The causal structure of the first order system (24-32) is 
simpler than expected: one can easily check from either 
(20,23) or (24-32) that 

dtiTi + 2Zi)= dt{2 Ei-Di-Ai + 2 Z^) = ... (33) 

(the quantities Ti here are just the space components of 
the quantities T ^ in eq. 20) so that one readily identi- 
fies three eigenfields propagating along the normal lines. 
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These "standing modes" are the only deviation of our sys- 
tem with respect to the wave equation pattern: all the 
remaining non-trivial eigenfields propagate along light 
cones. If we select a specific space direction, along a 
given unit vector Uk, then these light cone eigenfields are 

K,j ± UkX'^ij , (trK - 2 e) ± UkA''. (34) 

A straightforward calculation shows then that the first 
order system (24-27) is symmetric hyperbolic. The "sym- 
metrizer" can be easily identified starting from quadratic 
positive-definite "Energy" functions which are conserved 
up to lower order terms. One such "Energy estimates" is 
provided by 

E = K'^Kij + X^'^Xkij + {trK - 2 6)^ (35) 
+ A'^Ak + (r'^ + 2 Z''){Tk + 2 Zk). 

We are currently working with finite difference algo- 
rithms that can take advantage of the symmetric hypcr- 
bolicity of the system (24-32) to increase the robustness 
of numerical simulations. Our preliminary results indi- 
cate that explicit "Energy" expressions of the form (35) 
can be very useful to devise stable boundary conditions, 
along the lines sketched in refs [8,37]. 
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